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Observations by the PHENIX and STAR collaborations suggest that a strongly coupled quark-gluon plasma 
is produced in heavy-ion collisions at the Relativistic Heavy Ion Collider (RHIC). After a brief introduction 
to heavy-ion physics, measurements of heavy-quark production in heavy-ion collisions and the modification 
of heavy-quark spectra by the QGP are presented. Measurements of the total charm cross-section in several 
different collision systems confirm that cc pairs are produced through parton hard-scattering in the initial stages 
of the collisions. Non-photonic (proxies for heavy quarks) are suppressed by a factor of ~ 5 in central Au 
+ Au collisions relative to p -|- p collisions. This is larger than most current theoretical predictions and has led 
to a re-examination of heavy-quark energy loss in the medium. The relative contributions of c and b decays to 
the non-photonic spectrum have been predicted by perturbative QCD calculations; STAR measurements of 
azimuthal correlation functions of non-photonic and hadrons agree with these predictions. 



1. Introduction 



The main purpose of the PHENIX and STAR ex- 
periments at the Relativistic Heavy Ion Collider is to 
study the properties of strongly interacting matter at 
high temperatures. Figure [T] shows a schematic the- 
oretical phase diagram of nuclear matter for various 
temperatures and baryon chemical potentials P, 
At normal baryon chemical potentials and low tem- 
peratures, strongly interacting matter exists as nu- 
clei or a gas of interacting hadrons. At higher tem- 
peratures, however, the quarks may become decon- 
fined due to asymptotic freedom; the degrees of free- 
dom of the system are not hadrons, but individual 
quarks and gluons. This state of matter, called the 
quark-gluon plasma (QGP), is thought to have ex- 
isted in the first few microseconds following the Big 
Bang, and may also have been produced in high- 
energy nucleus-nucleus collisions such as those at 
RHIC at Brookhaven National Laboratory and the 
SPS at CERN. 



Quark-Gluon Plasma 
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Figure 1: Schematic phase diagram of strongly interacting 
matter. On the vertical axis is temperature and on the 
horizontal axis is baryon chemical potential, ^baryon- 



Lattice QCD calculations predict a sudden increase 
in the number of degrees of freedom of strongly inter- 
acting^ matter near a critical temperature Tc ~ 170 
MeV (l| . Whether or not there is a true phase transi- 
tion and its exact nature has yet to be determined. It 
has been proposed that a critical point may exist in 
the phase diagram, and that for higher baryon den- 
sities, a first-order phase transition may be observed 
[l|, 0]- A search for this critical point may be con- 
ducted in the next few years at RHIC by reducing the 
collision CM energy to lower values (^^SGeV) than 
the collider's present operating range (22 GeV to 200 
GcV). 

The collision of two nuclei introduces a large 
amount of energy into a region of space approximately 
the size of a nucleus for a short period of time. Lat- 
tice calculations predict a critical energy density of 
« 700 MeV/fm3 needed for QGP formation 0. The 
highest-energy RHIC collisions reach an energy den- 
sity of at least 4.9 GeV/fm^, seven times the criti- 
cal energy density [|. The QCD vacuum is "melted" 
and a quark-gluon plasma is produced. It is believed 
that the QGP quickly reaches thermal equilibrium, 
expands, and cools for a few fm/c until the transition 
temperature is reached. At this point, the partons be- 
come confined into hadrons. The hadron gas expands 
and the hadrons scatter inelastically until chemical 
freeze-out. The hadron gas continues to expand for 
a few more fm/c with elastic hadron- hadron interac- 
tions until thermal freeze-out, after which hadronic 
interactions are negligible Q . The resulting shower of 
particles can be detected in detector systems at RHIC. 

Many phenomena in heavy-ion physics depend upon 
the degree of overlap between the two colliding nu- 
clei, called the centrality of the collision. If the dis- 
tance between the centers of the nuclei (impact pa- 
rameter) is small, the overlap between the nuclei is 
large. Such collisions are called "central" events. A 
peripheral event has a small overlap and the impact 
parameter approaches the sum of the radii of the two 
nuclei. Npart is the number of nucleons that partici- 
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pate in collisions and Nunary is the number of binary 
collisions between those participating nucleons. Npart 
and Nfjinary are large for central events and small for 
peripheral events. The centrality of an event is es- 
timated using the multiplicity of charged tracks at 
mid-rapidity. Central collisions are characterized by 
a higher charged-particle multiplicity than peripheral 
events. Npart and Nunary are estimated using the 
Glauber model of nucleus- nucleus interactions 0, @| . 

The ratios of various particle yields have been fit 
with statistical thermal models, which indicate that 
at the time of chemical freeze-out, the system is ther- 
malized and has a temperature T = 170 — 180 MeV 
for RHIC collision energies A non-central colli- 
sion will have a spatially asymmetric overlap region, 
roughly elliptically shaped. The spatial asymmetry 
translates into a momentum-space asymmetry in the 
spectra of final state particles, with hadrons emitted 
preferentially near the reaction plane (the plane con- 
taining the beam axis and the impact parameter vec- 
tor). Measurements of this phenomenon (called el- 
liptic flow) are well described by ideal hydrodynam- 
ical models, indicating that the state of matter pro- 
duced in the collision has a very low viscosity and is 
a nearly perfect liquid @, i, [13. AU RHIC experi- 
ments have observed the production of back-to-back 
di-jets in p-\-p, d-\- Au, and peripheral Au -I- Au colli- 
sions. In central Au -f Au collisions, however, di-jets 
are not observed: the away-side jets appear to have 
been quenched by the medium pj] . [T^ . 
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Figure 2: Nuclear Modification Factor, Raa, for light- 
flavor hadrons and direct photons, tt'' and r) are sup- 
pressed by a factor of ~ 5 in central Au + Au collisions at 
y/SNN = 200 GeV relative to p + p collisions at the same 
energy. Direct photons are not suppressed. The curve was 
generated using a theoretical model of parton interactions 
with the QGP. The data are used to constrain the value 
of the gluon density [dNg/dy) in that model. 

Hard-scattering processes (those with large momen- 
tum transfer) between the incoming nucleons are ex- 
pected to account for 30-50% of the particles pro- 
duced in a nucleus-nucleus collision. Hard scatter- 
ing often produces particles with high transverse mo- 



mentum [pt > 2 GeV/c). Partons passing through 
a quark-gluon plasma are expected to lose energy to 
the medium through gluon radiation and collisions 
with other partons in the medium [l^ [l3| . The pro- 
duction of a QGP in a nucleus-nucleus collision will 
therefore cause a depletion of high-pr particles rela- 
tive to proton-proton collisions, which do not produce 
a medium. The standard measure of the effect of the 
medium on particle yields is the nuclear modification 
factor Raa- Raa is the ratio of a particle yield in 
nucleus- nucleus collisions to the yield produced in p+p 
collisions. The ratio is scaled by I / (Nbinary) , where 
(Nbinary) is the average number of nucleon-nucleon 
collisions in a nucleus-nucleus collision. 



Raa 



d^N{A+A) 
dprdy 



{N binary 



d?N{p+p) 
dpTdy 



(1) 



If no medium is produced, then a nucleus- nucleus col- 
lision can be viewed as an incoherent superposition 
of nucleon-nucleon collisions and Raa will be unity. 
Deviations from unity indicate the effects of nuclear 
matter and the quark-gluon plasma on particle yields. 
Central collisions produce a larger medium than pe- 
ripheral collisions, which should result in a greater 
suppression of high-pT particles. 

Figure [2] shows measurements of Raa for direct 
photons, tt", and rj by PHENIX. Raa is unity for 
direct photons [3], indicating that there is no sup- 
pression of direct photons in nucleus-nucleus colli- 
sions. This is expected since photons do not interact 
strongly and should not be affected by the presence 
of a QGP. However, the hadrons [H, [13) [ill are sup- 
pressed by about a factor of 5 at high pT- Also shown 
is a theoretical prediction of light-flavor-hadron sup- 
pression from the GLV model of parton-QGP inter- 
actions [igl . [20j . dNg/dy is the gluon density, a pa- 
rameter in the GLV model related to the opacity of 
the medium. The observed light-flavor-hadron sup- 
pression was used to constrain the value of this pa- 
rameter, giving dNg/dy « 1100. Raa has also been 
measured for heavy-flavor decay products and com- 
pared to models; this will be discussed in subsequent 
sections. These and other measurements indicate that 
the matter produced in RHIG collisions is a strongly 
coupled quark-gluon plasma (sQGP) [2l|. 



2. Experiments 

The experiments described in these proceedings 
are conducted at the Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven National Laboratory on Long 
Island, NY [22|. Ions are accelerated from the tan- 
dem Van de Graaff accelerators or the proton sources, 
through the AGS Booster and the AGS, and are in- 
jected into RHIC. Electrons are stripped off at sev- 
eral locations along the way. RHIG consists of two 
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synchrotron rings, 3.8 km in circumference. The two 
counter-circulating beams of ions intersect in each of 
the six interaction regions. The PHENIX and STAR 
detectors sit at two of these interaction regions. At 
RHIC, physicists can create and maintain beams of 
ions ranging from protons (both polarized and unpo- 
larized) to the heaviest nuclei. RHIC can be used to 
collide protons with a center-of-mass collision energy 
up to — 500 GeV and a luminosity of ~ 10^^ 
cm~^s~^. RHIC can also collide ions ranging in mass 
from deuterons to gold nuclei with CM energy per 
nucleon pair ^/sJ^N < 200 GeV and a luminosity of 
- 10^6 cm~2s-i [2|. 



Coils BEMC 




Figure 3: The STAR detector. Measurements described 
below use the Time Projection Chamber (TPC), Time-of- 
Flight detector (TOF), and the electromagnetic calorime- 
ter (EMC). 

Figure [3] shows a diagram of the STAR (Solenoidal 
Tracker At RHIC) detector (H]. The barrel and 
forward time-projection chambers (TPC and FTPC) 
record particle trajectories inside a room-temperature 
solenoidal 0.5-T magnet. The TPC covers pseudo- 
rapidity |?7| < 1.8, while the FTPCs cover 2.5 < 
\ri\ < 4. After passing through the TPC, particles 
enter the barrel electromagnetic calorimeter (BEMC) 
or the endcap electromagnetic calorimeter (EEMC, 
not shown). The STAR calorimeters together cover 
pseudorapidity -1 < 77 < 2. The Shower Maximum 
Detectors (SMDs), located approximately 5 radiation 
lengths inside each EMC tower module, provide ad- 
ditional particle identification based on the shape of 
the electromagnetic shower produced in the calorime- 
ters. The SMDs allow shower shapes to be measured 
to high precision (Ary = 0.007, Acj) = 0.007 rad). A 
silicon vertex tracker (SVT) covers |?7| < 1 between 
the beampipe and the TPC, providing accurate par- 
ticle tracking near the collision vertex. A prototype 




West Beam View 




Figure 4: The PHENIX detector in 2003. Measure- 
ments described below use the Pad Chambers (PC), Drift 
Chambers (DC), Time Expansion Chambers (TEC), Ring- 
Imaging Cherenkov detectors (RICH), Time-of-Flight de- 
tector (TOF), and the lead-scintillator (PbSc) and lead- 
glass (PbGl) electronmagnetic calorimeters. 



Time-of-Flight detector (TOF) has been installed out- 
side the TPC covering — 1 < < and A(j) = 0.04 rad. 
The TOF provides a precise measurement of particle 
velocity. Plans call for the TOF to be extended to full 
azimuthal coverage over pseudorapidity \r]\ < 1. A 
full description of the STAR detector is given in [2^ . 

Figure [4] shows a diagram of the PHENIX 
(Pioneering High-Energy Nuclear Interaction 
Experiment) detector. [2J| The two central spec- 
trometer arms sit in an axial magnetic field and 
cover I77I < 0.35 and tt/2 each in azimuth. Particle 
tracking is provided by three sets of pad chambers 
and (in the east arm) the time expansion chamber. 
The Ring-Imaging Cherenkov detectors and the 
Time-of-Flight detector provide particle identifica- 
tion. Beyond these detector systems sit lead-glass 
and lead-scintillator calorimeters. The two forward 
muon spectrometer arms sit in radial magnetic fields. 
They consist of drift chambers for precision tracking 
and muon identifiers. The muon identifiers are made 
up of alternating layers of steel absorber plates and 
streamer-tube tracking layers. A full description of 
the PHENIX detector, including the azimuthal and 
pseudorapidity coverage of each detector subsystem, 
is given in 24]. 
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3. Heavy Flavors 



Due to the high luminosity at RHIC, even parti- 
cles with comparatively low production cross-sections, 
such as c and h quarks, can be used to probe the 
strongly interacting matter produced. In nucleus- 
nucleus collisions, the dominant production mecha- 
nism for heavy quarks is gluon-gluon fusion in the ini- 
tial hard scattering of nucleons [23 . Thermal produc- 
tion of heavy-quark pairs in the medium is believed 
to be negligible after the initial stages of the colli- 
sion and the number of heavy quarks is essentially 
"frozen." Therefore, heavy quarks are probes sensi- 
tive to all stages in the evolution of the QGP, from its 
initial formation to hadronization and freeze-out. 

The presence of a QGP is expected to lead to 
the dissocation and suppression of heavy quarkonia 
through Debye screening of color charges. Some the- 
oretical calculations predict sequential dissociation of 
heavy quarkonia, with the more weakly bound reso- 
nances dissociating at lower temperatures [1^ . If this 
is true, measurements of heavy quarkonium suppres- 
sion could lead to a measurement of the initial tem- 
perature of the medium. PHENIX measurements of 
Raa for J/ V' 113, m, nil show less suppression than 
was expected based on SPS data. This may be due 
to the regeneration of J /ip through the recombination 
of dissociated charm (anti)quarks into new cc pairs 

H [Mill Sill. 

While all partons lose energy to the medium 
through gluon bremsstrahlung, the intensity and an- 
gular distribution of the gluon radiation is predicted 
to depend upon the mass of the radiating parton. In 
the small-angle approximation, the angular distribu- 
tion of gluon radiation from heavy quarks differs from 
the bremsstrahlung spectrum for massless partons by 
the pre- factor given in equation [2l 

dP (heavy quark) Mq 1 2 dP [light quark) 

(2) 

where Mq and Eq are the heavy-quark mass and 
energy, respectively jsH]. This suppression of gluon 
radiation from heavy quarks is strongest for angles 
less than Mq/Eq. The suppression of small-angle 
gluon radiation is called the dead cone effect. Be- 
cause of the dead cone effect, heavy quarks should 
lose less energy to the medium through gluon radia- 
tion than light quarks. Partons can also lose energy 
through collisions with other partons in the medium. 
For heavy quarks traversing the medium, the amount 
of energy lost through collisions may be comparable 
to the amount lost through gluon radiation (isl . IT3 | . 



4. Measurements of Open Heavy Flavors 

The RHIC experiments have studied several dif- 
ferent heavy- flavor decay channels, including the 
hadronic decays of mesons to pions and kaons, 
and the semileptonic decays of heavy-flavor hadrons 
to muons and electrons. 

4.1. Direct Reconstruction of Decays 

The STAR collaboration has found the yields of D° 
and D'^ mesons [s^ by reconstructing the D'^{D'^) 

+ TT^ decays, which have a branching ratio of 
3.83% [13] ■ The pions and kaons are identified by 
their energy loss in the Time Projection Chamber. 
STAR cannot reconstruct the full decay topology since 
ct{D°) = 124 ^m [33 and the TPC does not have suf- 
ficient track projection resolution to distinguish 
decay products from tracks coming directly from the 
primary collision vertex. The D'^ invariant mass spec- 
trum was obtained by pairing each kaon with oppo- 
sitely charged pions from the same event. The com- 
binatorial background was estimated through event 
mixing techniques and subtracted. Figure [5^ shows 
the Ktt invariant mass spectrum for lyl < 1 in Au -I- 
Au collisions at ^^snn = 200 GeV [3^. A clear 
peak is visible. Because of the small branching ratio 
for this decay and the lack of a dedicated trigger, this 
analysis is limited by statistics to pt < 3 GeV/c. 

4.2. Decays to Muons 

The STAR collaboration has identified low-px 
muons (0.17 GeV/c < pr < 0.25 GeV/c) through 
measurements of energy loss in the Time Projection 
Chamber and = (p/iPj))"^ in the Time-of-Flight 
detector and the TPC. Figure [5)d shows the dis- 
tribution of particles after energy- loss selection [soj . 
The muon and pion peaks are clearly visible. In ad- 
dition to muons from heavy flavor decays, the muon 
peak contains a large number of muons from pion and 
kaon decays; it was necessary to remove these "back- 
ground" muons. HIJING [40] was used to simulate the 
DCA (distance of closest approach) of muon tracks 
with the primary collision vertex (see Figure [5j;). The 
DCA distribution for muons from charm decays has 
a maximum much closer to zero than the distribution 
for muons from pion and kaon decays. The observed 
DCA distribution was fit with a linear combination of 
these two simulated DCA distributions to obtain the 
contribution of charm-decay muons to the total muon 
yield. [39] 

4.3. Decays to Single Electrons 

Both the PHENIX and STAR collaborations have 
also studied the spectrum of single electrons, i.e. 
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Figure 5: STAR heavy-flavor identiflcation methods, (a) D'^ meson reconstruction in 200-GeV Au+Au collisions. The 
Ktt invariant mass distribution is shown after the subtraction of the combinitorial background (estimated using event- 
mixing). The residual background is approximated by a linear function and subtracted, (b) distribution measured 
by TOF detector in 200-GeV An + Au collisions. The particles have passed the energy-loss cut used to identify low-pr 
muons. (c) DCA distributions used to identify charm-decay contribution to the inclusive muon spectrum (see explanation 
in text). 



produced with (anti) neutrinos in weak decays. The 
spectrum of single is expected to be dominated 
by heavy-flavor semileptonic decays (e.g. D'^ 
e+ + + Vq). The main sources of background to 
the single signal are e+e~ pairs from tt*^ and rj 
Dalitz decays and photon conversions [sj- For this 
reason, single are called non-photonic e^, while 
the background are called photonic e"^. The de- 
cays of vector mesons (e.g. p, cj), and to) make small 
contributions to the photonic background. Back- 
ground from photon conversion is less significant at 
PHENIX due to the reduced amount of material rel- 
ative to STAR. Ke3 decays [K^ ^ n° + + Ue{ve) 
and if£ — > + + i/eii^e)) make small contributions 
to the single signal. 

The STAR collaboration identifies using two dif- 
ferent methods. For pT < 3.5 GeV/c, are identifled 
through measurements of energy loss and with the 
Time Projection Chamber and Time-of-Flight detec- 
tor (similar to the method of muon identification de- 
scribed above) [11]. For pr > 1-5 GeV/c, e"^ are iden- 
tified through measurements of TPC energy loss, the 
energy E deposited in the electromagnetic calorime- 
ter, and the shape of the electromagnetic shower mea- 
sured in the shower maximum detector [4lj. The use 
of a trigger allows this measurement to extend up to 
Pt ~ 8 GeV/c. An energy- loss cut of approximately 
3.5 keV/cm < dE/dx < 5 keV/cm is used to identify 



(the exact bounds are varied slightly depending on 
track momentum and event charged-track multiplic- 
ity). Compared to hadrons, produce larger show- 
ers and deposit more of their energy in the EMC. The 
ratio p/ E has a maximum near 1 for (see Figure 
[S] 42]) and additional identification is provided by 
cuts on the shower size in the SMD. The remaining 
hadron contamination is « 2% at ~ 2 GeV/c and 
« 20% at pT 8 GcV/c. STAR identifies photonic 
(background) through invariant-mass reconstruc- 
tion of e+e^ pairs. The invariant-mass distribution 
of e+e~ pairs from photon conversions and 7r° and rj 
Dalitz decays has a maximum near mass. When 
a cut of Minv{e^ e~) < 150 MeV/c^ is used, photonic 
electrons are identified with an efficiency around 70%. 

The PHENIX collaboration identifies using the 
Ring-Imaging Cherenkov detector, measurements of 
the shower shape in the electromagnetic calorimeter, 
and a cut on the energy-to-momentum ratio [4^ . . 
Photonic are identified using two different meth- 
ods. In the "cocktail subtraction" method, the spectra 
of from various sources of background are simu- 
lated. Measured yields of 7r°, ij, direct photons, and 
other sources of background are used as input for the 
simulation generator. In the "converter subtraction" 
method, a photon converter (a thin brass sheet of 
1.67% Xq) is inserted around the beam pipe. ANg, 
the increase in the yield due to the converter, is 
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Figure 6; p/E distributions for and hadrons in 200- 
GeV p+p collisions at STAR, and hadrons are identified 
through their energy loss in the Time Projection Cham- 
ber. The hadron distribution has been scaled to match 
the distribution at high values of p/E. This provides 
an estimate of the hadron background that remains in the 
peak. 



measured. GEANT simulations are used to determine 
i?-),, the fractional increase in the photonic yield 
caused by the converter {R^ « 2.3 in the 2006 p + p 
run). Knowledge of ATVg and R^ allows the photonic 
yield to be determined and removed from the in- 



clusive e yield, 
with no converter present, 



Where is the photonic 



yield 



(3) 



The non-photonic yields measured using the cock- 
tail and converter subtraction methods are consistent 
with each other. 



shown [45]. The cross-section is divided by {Nbinary), 
the average number of binary collisions for the given 
collision system. Within each experiment, the charm 
cross-section scales with the number of binary colli- 
sions. This is a confirmation that charm is indeed pro- 
duced through initial parton hard-scattering and that 
charm production through other mechanisms (such as 
thermal production in the QGP) is not significant. 
Note that the cross-sections measured by STAR are 
higher than those measured by PHENIX by a factor 
of ^2. This is still under investigation. 

Due to the large quark masses, the hard-scattering 
processes that produce heavy flavor can be calculated 
using perturbative QCD [4^. The most advanced 
perturbative calculation scheme is the Fixed-Order 
plus Next-to-Leading-Log-resummed approximation, 
or FONLL. A FONLL prediction for the charm cross- 
section [4^, im is shown as prediction (I) in Figure 
[5^. The PHENIX data are consistent with this pre- 
diction; the STAR data are greater than the FONLL 
prediction by a factor of ^5 and sit well above the 
upper uncertainty bound. However, the charm cross- 
section predictions are sensitive to the number of ac- 
tive flavors, the choice of scale, and the parton densi- 
ties in the collision system. A new calculation by R. 
Vogt [ill (prediction (II) in Figure [SK) indicates that 
the uncertainties on the perturbative calculation may 
be larger than previously thought. The total charm 
cross-section measured by STAR is consistent with the 
new perturbative calculation. 

Figure shows the ratios of the measured non- 
photonic e± yields [H, IHl, ^ to the FONLL predic- 
tion as functions olpT for p-\-p collisions at — 200 
GeV. The dashed horizontal lines indicate the ratio 
of the measured cross-sections to the FONLL cross- 
section. Figure [8)3 indicates that STAR's disagree- 
ment with PHENIX (by a factor of -2) and FONLL 
(by a factor of ^--^5) exists even at high pT- However, 
the FONLL prediction does describe the shape of the 
STAR and PHENIX non-photonic spectra well. 



5. Total Charm Cross-Section 

STAR determines the total charm cross-section, CTce, 
for each collision system through a combined fit of 
the -D", muon, and non-photonic measurements 
described above (e.g . Figure [7^ for 200-GeV Au 
+ Au collisions) [3^. PHENIX determines the to- 
tal charm cross-section from the measurement of the 
non-photonic yield (e.g. Figure [Tls for 200-GeV 
p -I- p collisions) 43]. Figure [5^ shows measurements 
of the scaled total charm cross-section by STAR and 
PHENIX ioT p + p d + An [3^, and Au -f 

Au [41], [13] collisions (in different centrality bins) at 
y^SNN = 200 GeV. A preliminary STAR measure- 
ment for Cu + Cu using only Z)*^ reconstruction is also 



6. Medium Modification of Non-photonic 
Spectra 

The cross-section discrepancy between PHENIX 
and STAR cancels in the nuclear modification fac- 
tor Raa, which is the scaled ratio of particle yields 
(see definition above). Figure [9^ shows Raa of non- 
photonic for Au -I- Au colhsions at ys n n — 200 
GeV for both PHENIX [H] and STAR ^,\^. Raa 
is plotted as a function of (Npart), the average num- 
ber of nucleons participating in collisions in a given 
centrality bin. The most central collisions have the 
highest values of {Npart)- The PHENIX and STAR 
Raa are consistent with each other across the range 
of {Npart) shown. Plotted in shaded bands are the val- 
ues of Raa measured by PHENIX for 7r° [l6( and by 
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Figure 7: Measurements used in calculation of total charm cross-section, (a) STAR 73", muon, and spectra for 200- 
GeV Au + Au collisions in various centrality bins, (b) PHENIX non-photonic spectrum for 200-GeV p + p collisions. 
Also shown are pQCD predictions (FONLL) for the c- and 6-decay contributions to the spectrum (see discussion below 
and Figure [TD]). 
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Figure 8; Comparisons of RHIC measurements to FONLL predictions, (a) Total charm cross-section (scaled by 
l/{Nbinary)) in different collision systems at y^sjvjv ~ 200 GeV. Theoretical predictions are also shown, (b) Ratios 
of non-photonic e"^ yields in central 200-GeV Au + Au coUisions for STAR and PHENIX to FONLL prediction (I) in 
central 200-GeV Au + Au collisions. STAR have been selected using measurements of energy-loss (dE/dx) in the 
TPC or in the TOF. 



STAR for charged light-flavor hadrons [13] ■ The mea- 
sured suppression of non-photonic is similar to the 
suppression observed for light-flavor hadrons. This 
was unexpected, as heavy quarks were expected to 
lose less energy in the medium than light quarks. As 
a result, non-photonic (proxies for heavy quarks) 
would be suppressed less than light-flavor hadrons. 



Figure ^jp shows a comparison of the non-photonic 
suppression measured by STAR and PHENIX to 
several theoretical models of heavy-quark interactions 
with a quark-gluon plasma. The data are from cen- 
tral Au + Au collisions at y^sjvjv — 200 GeV. In 
d + Au collisions at this energy a ~ 50% enhance- 
ment of non-photonic is observed (data [4l| not 
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Figure 9: Nuclear modification factor {Raa) for non-photonic in Au + Au collisions at 200 GeV. (a) Raa vs. (Npart)- 
Also shown are light-flavor hadron Raa data, (b) Raa vs. pr for central collisions. The data are compared to theoretical 
predictions. 



shown). The enhancement is explained by the Cronin 
effect (multiple scattering in normal nuclear matter). 
Therefore, the strong non-phonic-e^ suppression ob- 
served in central Au -I- Au collisions does not appear 
to be due to such cold nuclear matter effects. The 
calculations that produce curves (I) and (II) include 
only energy loss through gluon radiation; these calcu- 
lations under-predict the suppression of non-photonic 
e± [5ll,[53|. dNg/dy is the gluon density in the DGLV 
model (curve (I)) and q is the time-averaged transport 
coefficient in the BDMPS model (curve (II)). These 
parameters are related to the opacity of the QGP, 
with higher values corresponding to a more opaque 
medium. Light-hadron suppression data are used to 
constrain the values of these parameters. The DGLV 
model with dNg/dy « 1000 and the BDMPS model 
with g « 10 GeV^/fm describe the observed light- 
hadron suppression well. 

In addition to energy loss through gluon radiation, 
partons can lose energy through collisions with other 
partons in the medium. Calculations [Tsj indicate 
that, for heavy quarks, collisional energy loss is as 
important as radiative energy loss. Curve (III) is gen- 
erated using the DGLV model (which produced curve 
(I)) including both radiative and collisional energy 
loss [1^. Curve (IV) is generated using the model 
of van Hees et al. '54], which includes heavy-quark 
energy loss through elastic coUsions in the medium 
and the formation of resonant D- and i3-meson states 
through quark coalescence. Curves (III) and (IV) 
seem to describe the data better than curves (I) and 



(II), which include only radiative energy loss, but still 
tend to under-predict the observed suppression of non- 
photonic e^, especially at high px- The model of Adil 
and Vitev [55| (curve (VI)) uses the fragmentation 
of heavy quarks and the subsequent dissociation of 
D- and i3-meson states in the medium to explain the 
suppression pattern of non-photonic e^. This model 
does describe the observed suppression well. 

Curves (I)- (IV) and (VI) were generated using 
FONLL predictions [i^l for the relative contributions 
of c- and b- quark decays to the non-photonic 
spectrum. Curve (V) is generated with the assump- 
tion that only c-quark decays contribute to the non- 
photonic spectrum and that the 6-quark contribu- 
tion is insignificant for the px range shown ^SSj . Curve 
(V) (which includes both radiative and collisional en- 
ergy loss) is consistent with the measured suppression. 
It is therefore important to disentangle the c and b 
contributions and determine their relative strengths; 
this will be discussed in the next section. 



7. Electron-Hadron Azimuthal 
Correlations 

Figure [10] shows FONLL calculations of the D- 
and D-meson-decay contributions to the non-photonic 
electron cross-section as functions of pr • The up- 
per and lower limits result from the variation of model 
parameters including quark masses and renormaliza- 
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tion scales. The charm contribution is dominant at 
low px- The bottom contribution becomes larger than 
the charm contribution around px ^ 5 GeV/c. 




5 10 15 



Pt (GeV) 

Figure 10: FONLL predictions for the contributions of D- 
and _B-meson decays to the mid-rapidity non-photonic 
cross-section at -/s — 200 GeV. 
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Figure 11: A(j> for pairs of non-pliotonic and hadrons. 
5.5 GeV/c < pT(e±) < 6.5 GeV/c andpT(ft) > 0.3 GeV/c. 
On the vertical axis is the number of charged tracks per 
trigger e*. 

To disentangle the charm and bottom contributions, 
the STAR collaboration has studied azimuthal corre- 
lations between non-photonic and hadrons. Fig- 
ure [Tl] shows the difference in azimuthal angle (A0) 
between non-photonic and hadrons in 200-GeV 
proton- proton collisions for one pt bin [56i] . Due to 
decay kinematics, there will be an azimuthal corre- 
lation between the leptons and hadrons produced in 
heavy-flavor semi-leptonic decays. Due to the larger 
B-meson mass, a S-meson can give more kinetic en- 
ergy to its decay products than a Z)-meson, resulting 
in a broader near-side (A0 ~ 0) correlation peak. The 
expected charm and bottom contributions to this dis- 
tribution are simulated with PYTHIA. Varying the 



■ Data points with fit range: -1.5 to 1.5 




2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 

(GeV/c) 

Figure 12: Fractional contribution of _B-meson decays to 
the non-photonic spectrum at y/s — 200 GeV. 

quark fragmentation functions inside PYTHIA does 
not significantly alter the shapes of the simulated 
distributions. The measured distribution is then fit 
with a linear combination of the simulated charm and 
bottom distributions to find their relative strengths. 
B / {B + D) , the fraction of the total non-photonic 
cross-section due to B-meson decays, is plotted in Fig- 
ure [T^l The measurement of B/{B + D) through e- 
h correlations ^561 is consistent with FONLL predic- 
tions. This is an indication that 6-quark suppression 
should be taken into account in describing the sup- 
pression of non-photonic at moderate px (see the 
discussion in the previous section). 



8. Conclusions 

Measurements of jet quenching and elliptic flow sug- 
gest that the state of matter created in a heavy-ion 
collision at RHIC ener gies is a strongly coupled quark- 
gluon plasma (sQGP) [2l| and is a near-perfect liquid 
PHENIX and STAR have calculated the total 
charm cross-section acs for p + p, d + Au, Cu -I- 
Cu, and Au -I- Au coUisions at ^/snn — 200 GeV 
[Si, SH, Si, El Ei]. The cross-section scales with 
the number of binary nucleon-nucleon collisions, an 
indication that charm quarks are indeed produced 
through initial parton hard-scattering. PHENIX and 
STAR have measured the nuclear modiflcation fac- 
tor Raa for light-flavor hadrons [11] and for non- 
photonic e*, which come predominantly from heavy- 
flavor semileptonic decays [4ll [isl . . The suppres- 
sion of non-photonic in central Au -I- Au collisions 
is greater than expected. It has been difficult for the- 
oretical models [5l|, [H, [Hi, [H, [5^1 to describe the 
suppression of light-flavor hadrons and non-photonic 
simultaneously. Perturbative QCD calculations 
(FONLL) [47I predict that the 6-decay contribution 
to the non-photonic becomes comparable to the 
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c-decay contribution sA, pr ^ ^ GeV/c. STAR mea- 
sures the difference in azimuth angle (j) for pairs of 
non-photonic and hadrons (56j . This distribution 
is fit with PYTHIA simulations of the expected c- 
and 6-decay contributions to determine their relative 
strengths; these data are consistent with the FONLL 
predictions. 
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